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Abstract 

For a system of coupled anharmonic oscillators we compare the conver- 
gence rate of the variational collocation approach presented recently 
by Amore and Fernandez (2010 Phys.Scr.81 045011) with the one ob- 
tained using the optimized Rayleigh-Ritz (RR) method. The mono- 
tonic convergence of the RR method allows us to obtain more accurate 
results at a lower computational cost. 

1 Introduction 

Since the 1930 s, the hnear Rayleigh-Ritz (RR) variational method [1] has 
been the primary tool for determining the low energy part of the spectrum 
of quantum systems. In this approach, the trial function is expanded into 
a finite set of basis functions, and the coefficients of the expansion are de- 
termined by the variational principle. The approximations provided by the 
eigenvalues of the RR matrix are upper bounds that approach the exact ener- 
gies monotonically from above on increasing the number of included functions 
of the basis [Tj. However, the matrix elements are given by integrals and the 
number of the basis functions required to obtain accurate results grows expo- 
nentially with the number of degrees of freedom, which renders the method 
computationally prohibitive for complex systems. There are basically two 
strategies to reduce the computational cost of the method: choosing the ba- 
sis functions in a sophisticated way or using approximations in calculating 
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the Hamiltonian matrix elements. Among various implementations of the 
first strategy the introduction of nonlinear parameters into the functions of 
the basis is one of the most effective The adjustable parameters may 
be fixed so as to improve the convergence for a particular state [3] or a 
group of states [4J. An example of the second strategy is the collocation 
method where the matrix elements are approximated by A^— point quadra- 
tures. The recently published paper on "Variational collocation for systems 
of coupled anharmonic oscillators" [5] combines both strategies mentioned 
above, demonstrating that by introducing adjustable parameters the con- 
vergence of the collocation method may be improved. However, it should be 
remembered that in contrast to the RR method, the collocation method does 
not have a variational character. Discretization in the collocation approach 
introduces approximation in the Hamiltonian matrix, as a consequence of 
which the results do not represent upper bounds to the exact bound-state 
energies. Introducing the adjustable parameters does not change this fact; 
therefore the method proposed by Amore and Fernandez should be termed 
optimized rather than variational. 



2 The model 

To illustrate the similarities and differences between optimization in the RR 
and collocation methods, we consider the example of the PuUen-Edmonds 
Hamiltonian 

that has been studied in Ref. |i5j. The parameter A represents the coupling 
between the two degrees of freedom. Amore and Fernandez observed that 
introducing a rotation parameter, the convergence of the collocation method 
may be improved further. The optimal rotation angle has been found to be 
close to 7r/4, which corresponds to the Hamiltonian 

Hamiltonians ([1]) and (|2]) have the same spectrum but the speed of conver- 
gence of the collocation method turns to be considerably higher in case of 
the application to the rotated one ([2]). 
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Below we consider the application of the RR method to both Hamilto- 
nians using basis functions (j)'^^{x,y) = (pni^)4>miy) constructed as tensor 
products of the functions (pni^) from a basis in one-dimensional (ID) space 
that depends on an arbitrary parameter a. The Mth order approximations 
are obtained with the use of M even-parity ID functions, which corresponds 
to diagonalization of the RR matrix of dimension M^. Optimization is per- 
formed by choosing the value of a so as to make the trace of that matrix 
stationary at each order of approximation [4j. In our procedure, we take the 
advantage of symbolic calculation in the Mathematica package by performing 
the numerical algebra exactly in terms of the ratio of integers, which allows 
determination of the eigenvalues to arbitrary accuracy. 

3 Trigonometric basis 

First we calculate the accurate values of ground-state energy of the coupled 
oscillators through the optimized RR method, using the basis of trigonomet- 
ric functions 



which are eigenfunctions of the infinite potential well of a half-width L that is 
taken as a nonlinear parameter. The same basis has been used by Amore and 
Fernandez in their collocation approach [5] . In table [H we compare the ap- 
proximations to the ground-state energy £o{l) and £o{2) for the Hamiltonian 
(II]) and ([2]), respectively, for different values of M. The optimal value of the 
nonlinear parameter Lopt that renders the trace of the given order RR matrix 
stationary is also shown. The underlined digits are those that coincide with 
the exact result. We observe that both the series of approximations converge 
monotonically from above to the exact value with increasing M, but the con- 
vergence for the Hamiltonian ([2]) is much quicker. This is similarly true at 
other couplings as can be seen in table [21 where the results obtained with 
M = 35 are compared for different values of A. We conclude that the same 
rotation that speeds up the convergence of the collocation method improves 
also the convergence of the RR method. 




(3) 
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M 


Loptil) 


eo(l) 


Lopti'^) 


^o(2) 


10 


2.627 


3.01970464 


3.015 


3.01917772606335945122 


15 


2.988 


3.01921402 


3.443 


3.01917771478547514131 


20 


3.278 


3.01918109 


3.787 


3.01917771477200287135 


25 


3.525 


3.01917810 


4.077 


3.01917771477196754393 


30 


3.741 


3.01917777 


4.332 


3.01917771477196738793 


35 


3.935 


3.01917772 


4.560 


3.01917771477196738692 



Table 1: RR approximations to the ground-state energy of Hamiltonians ([T]) 
and (|2]) with A = 10 obtained using M functions from the trigonometric 
basis with optimized quarter period Lopt- The under hned digits are those 
that coincide with the exact result. 



A 


Loptil) 


eo(l) 


Lopti'^) 


eo(2) 


5 


4.41 


2.653909778 


5.10 


2.65390977795321535349058 


10 


3.93 


3.019177723 


4.56 


3.01917771477196738692058 


100 


2.69 


5.460971628 


3.12 


5.46097039792335813076119 


1000 


3.66 


11.23244729 


2.13 


11.2324392672100309737449 


10000 


1.83 


23.94601185 


1.47 


23.9459896278197863410786 



Table 2: RR approximations to the ground state-energy of Hamiltonians ([T]) 
and ([2]) for several values of A obtained using M = 35 functions from the 
trigonometric basis with optimized quarter period Lopt- 



4 Harmonic oscillator basis 

For comparison, we applied the optimized RR method using a different basis, 
namely the set of eigenfunctions of the harmonic oscillator with an arbitrary 
frequency Vt that are given by 



-)'''^if.(v^a:)e"¥. (4) 



The approximations to ground-state energy obtained with the harmonic os- 
cillator basis are presented in tables |3] and |H Comparison with the results 
obtained with the trigonometric basis in tables [1] and |2] shows that the con- 
vergence is much faster for the harmonic oscillator eigenfunctions, which may 
be attributed to the fact that they fulfill the proper boundary conditions at 
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infinity. The convergence for tlie rotated Hamiltonian is mucli better 
tlian for the original one ([T]) and the obtained results agree with the best 
available calculations of the perturbative-variational method performed to 
an accuracy of 15 digits For Hamiltonian ([2]) with A = 5 the precision of 
our calculations with M = 35 harmonic oscillator functions reaches 30 digits. 
For larger values of A, a similarly high precision can be easily obtained with 
reasonably increased M, which would be not an easy task for the original 
Hamiltonian ([1]). 



M 


^^opt(l) 


^o(l) 


^^opt(2) 


^o(2) 


5 


3.65 


M2083109727 


3.04 


3.0192115512360457661267360 


10 


4.64 


3.01920677896 


3.77 


3.0191777151455247319696151 


15 


5.32 


3.01917868488 


4.30 


3.0191777147719932650455193 


20 


5.86 


3.01917776226 


4.71 


3.0191777147719673955033821 


25 


6.31 


3.01917771779 


5.07 


3.01917771477196738691 95328 


30 


6.71 


3.01917771501 


5.38 


3.0191777147719673869116933 


35 


7.06 


3.01917771479 


5.65 


3.0191777147719673869116789 



Table 3: RR approximations to the ground-state energy of Hamiltonians ([T]) 
and ([2]) with A = 10 obtained using M functions from the harmonic oscillator 
basis with optimized frequency flopt- 



A 


^^opt(l) 


£o(l) 


^opt(2) 


eo(2) 


5 


5.63 


2.6539097779536 


4.51 


2.6539097779532153534905698061 7 


10 


7.06 


3.0191777147930 


5.65 


3.01917771477196738691167893635 


100 


15.13 


5.4609704165739 


12.08 


5.46097039792335524 182772613188 


1000 


32.56 


11.232439458359 


25.98 


11.2324392672098516856244029369 


10000 


70.14 


23.945990215098 


55.94 


23.9459896278189396640103254023 



Table 4: RR approximations to the ground state-energy of Hamiltonians ([I]) 
and ([2]) for several values of A obtained using M = 35 functions from the 
harmonic oscillator basis with optimized frequency Qopt- 
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5 Comparison of the results and conclusion 



In this paper, we studied the convergence of the optimized RR method for 
the ground-state energy of the coupled anharmonic oscillators using two ba- 
sis sets: trigonometric functions and harmonic oscillator eigenfunctions. The 
comparison of the effectiveness of the methods with that of the optimized 
collocation method fSj is summarized in figure [1] by plotting the precision 
achieved in the computation performed with Mathematica ver. 7.0. The 
precision calculated as a difference between the approximate and exact re- 
sult 5eo is shown in a log-log plot as a function of CPU time. An important 
observation is that the rotation that has been shown to improve the con- 
vergence of the collocation method ^ accelerates also the convergence of 
both variants of the RR method. It is seen that far greater accuracy in the 
ground-state energy is obtained by using the rotated Hamiltonian (|2]) than by 
increasing the size of the basis. The optimized RR method is more effective 
with the harmonic oscillator basis than with the trigonometric one, but in 
both cases the convergence is quicker than that of the optimized collocation 
approach. This is opposed to the common understanding that the computa- 
tional cost of the collocation method is less then that of the RR method [?]. 
Our results show that the RR method with an appropriately optimized ba- 
sis offers a similar accuracy to the optimized collocation method with much 
lower computational cost. 

Another important advantage of the RR method is its variational char- 
acter, which allows us to state that the lower the achieved result, the closer 
it is to the exact one. Discretization in the collocation approach facilitates 
calculation of the Hamiltonian matrix but destroys the variational behavior. 
Normally the monotonic convergence from above is recovered when the num- 
ber of grid point is sufficiently large but this is not automatically granted. 
For instance we observe that the approximations obtained by the collocation 
method from the rotated Hamiltonian ([2]) in table 4 of [5] are below the 
accurate ground-state energies ( [6] and our table H]) and they decrease with 
increasing M, which suggests that they do not converge to the proper result. 
This is a warning that when using the collocation method one has to be very 
cautious about deriving accurate results. 

We believe that our conclusions are relevant not only for the basis sets 
considered in this work but also for other types of basis that are used in 
higher-dimensional problems. 
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Figure 1: Precision achieved by the optimized collocation (COLL) and RR 
method with trigonometric (TRIG) and harmonic oscillator (HO) eigenfunc- 
tions for the Hamiltonians ([T]) and ([2]) with A = 10 as functions of CPU 
time. 
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